from odl.discr.grid import sparse_meshgrid
import numpy as np

scale = 1.0
N = 32
t = np.linspace(-scale, scale, N)
grids = sparse_meshgrid(t, t, t)
for i, grid in enumerate(grids):
    print(i, grid.shape, grid)
# cube = grids[0] + grids[1] + grids[2]
# print('cube', cube.shape, type(cube))
